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Abstract 



The vortex state at the lower critical field, H c i, in clean anisotropic su- 
perconductors placed in an external field tilted with respect to the axis of 
anisotropy (c-axis) is considered assuming two possible arrangements: dilute 
vortex- lines or dilute vortexdine chains. By minimizing the Gibbs free en- 
ergies in the London limit for each possibility we obtain the corresponding 
lower critical fields as a function of the tilt angle. The equilibrium config- 
uration at H c i for a given tilt angle is identified with that for which H c \ is 
the smallest. We report results for parameter values typical of strong and 
moderate anisotropy. We find that for strong anisotropy vortex-line chains 
are favored for small tilt angles (< 7.9°) and that at 7.9° there is coexistence 
between this configuration and a vortex-line one. For moderate anisotropy we 
find that there is little difference between the vortex-line and the vortex-chain 
lower critical fields. 
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Motivated by the discovery of high-T c superconductivity in the cuprates, calculations 
of vortex properties in uniaxially anisotropic superconductors were carried out by several 
workers [[j]|§. One surprising result is that the interaction between a pair of straight vortex 
lines parallel to each other, tilted and coplanar with respect to the c-axis is attractive at 
distances larger than a value of the order of the penetration depth, as was first found by 
Grishin, Martynovich and Yampol'skii As a consequence, a chain of such lines has 

lower energy than the same vortex lines placed far apart. On the basis of this result it was 
suggested that at if cl the mixed state consists of a vanishing small density of vortex-line 
chains instead of vortex lines 

However, this is not correct in general for tilted vortex lines. The reason is that in order 
to determine the H c \ configuration it is necessary to minimize the Gibbs free energy in the 
limit of vanishing vortex density for each candidate vortex arrangement and, from it, obtain 
the corresponding if cl . The equilibrium configuration is the one with the smallest H c \. In 
this paper we carry out this calculation in detail assuming two possible arrangements: dilute 
vortex lines (DVL) or dilute vortex-line chains (DVLC). 

The Gibbs free-energy minimization for the DVL was carried out by Sudb0, Brandt and 
Huse [0] in the London limit with the goal of studying coexistence of vortex states at H c \. 
They find that, over a range of parameter values, there is one external field tilt angle for 
which two DVL states, differing by the vortex-line orientation with respect to the c-axis, 
coexist at H c i. In this paper we also investigate how these results are modified by vortex-line 
chain states. 

To carry out the Gibbs free-energy minimization we proceed along lines similar to those 
developed by Sudb0, Brandt and Huse J7j. We obtain H c \ as a function of the external field 
tilt angle for parameter values typical of strong and moderate anisotropy, both for DVL 
and DVLC, and determine which one of these is the equilibrium one at H&. For strong 
anisotropy we predict coexistence between DVL and DVLC at a particular external field tilt 
angle. 

We consider a bulk uniaxially anisotropic superconductor placed in a magnetic field H, 
with magnitude H and tilted with respect to the c-axis by an angle a. We assume that the 
vortex lines are straight and point in a direction making an angle 9 with the c-axis. 

In the London limit the superconductor is characterized by the penetration depths 
and A c , for currents parallel to the ab-plane and to the c-direction, respectively. The free 
energy per unit length of a generic arrangement of these vortex lines, F, can be written as 
the sum of pairwise interactions 

^=^E/(r*-r,), (1) 

where rj is the i-th line position vector in the plane perpendicular to the vortex line direction 
and /(r) is the Fourier transform of 

1 4- A 2 P 

f(V\ = P" 2 3(k) 1 + A e K (o\ 

J[ ) (l + A^Xl + A^ + A^) ' U 

where A 2 , = A 2 b sin 2 9 + A 2 , cos 2 9, x and y refer to two directions perpendicular to one another 
and to the vortex lines, with x coplanar with the c-axis, and g(k) is the vortex core cutoff 
function g,§||. 
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For both DVL and DVLC, F can be cast in the form 



F (N v ,9)=n v Ae(9) 



(3) 



where n v is the vortex-line density and A is the sample area perpendicular to the vortex 
lines direction. For the DVL we can neglect interactions between lines, so that s{9) = e s f{6), 
the vortex-line self-energy (e s f(8) = ($o/87r)/(0)). For DVLC we can neglect interchain 
interactions, so that e{9) = £ s f(9) + £ c h(Q), where £ c h(6) is the energy of interaction per 
vortex-line of a chain running along the x-direction (e c h(6) = ($q/87t) J2 n ^o f( x = n Li,y = 
0), L\ being the vortex-chain period). 

In order to calculate H c i(a) we have to minimize the Gibbs free energy with H and 
the sample volume fixed. The Gibbs free energy per volume for the above described vortex 
states is given by 



G (B,9;H,a) = ^-[e{9) - ^Hcos (9 - a)] 

<3>0 47T 



(4) 



where B = &on v is the modulus of the magnetic induction. 

We adopt for e s f(9) the expression derived by Sudb0 and Brandt |I0|| , using an elliptic 
core cutoff function, 



/m V w \ A^cos 2 ^ 7 2 (A^ + A^). 
eAO) = e -[H K1 ) + os2Q+ In 



(5) 



where Eq = ($o/47rA a fe) 2 , 7 = A a b/A c , k = \ a b/£,ab, £,ab being the ab-plane coherence length. 
This line energy was derived from the anisotropic Ginzburg-Landau theory using the Klemm- 
Clem transformations |ll| , [l2| . The elliptical core cut off has semimajor axis and semimi- 
nor axis cos 2 9 + £ 2 sin 2 9)~ 1 / 2 . 

Taking the derivatives of Go with respect to B and 9 and equating to zero we obtain, 
respectively 



e{9) - ^Hcos(9-a) = 



(6) 



!!g> + £™, ( ,- a) = o. (7, 

The first equation is the familiar condition Go = at H ci , whereas the second results from 
torque balance. 

To solve this system of equations we first obtain H from Eq. (^) and substitute it in Eq. 
(0). The following relationship between a and 9 results 

Finally, we obtain H c \(a) by substituting the function 9(a) obtained by inverting Eq. (||) 
into Eq. (§) 
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Hel{a) = *zm . ( 9) 

$ cos [9(a) - a] 

As we shall see shortly, 9(a) may be a multivalued function. Multivalued solutions for if cl in 
anisotropic materials have been first noticed by Klemm and Clem ||1 1|| . If this is so, according 
to Eq. (^), for each a there are several, generally distinct, H c \(a) and corresponding vortex 
configurations satisfying the minimization conditions. The physical solution is the one with 
the smallest H c i(a). 

Now we obtain numerically the solutions of these equations for the following choices of 
parameters, typical of strong and moderate anisotropy: 
Moderate anisotropy: k — 50, 7 = 1/5 
Strong anisotropy: K = 10, 7 = 1/^/200 
The results are as follows. 

i) Dilute arrangement of vortex lines (DVL). As expected, our results are identical to 
those obtained by Sudb0, Brandt and Huse For strong anisotropy there is a region where 
9(a) is multivalued: for each a there are three values of 9 and of H c ±(a) as shown in Fig. |]. 
The physical solution for each a is the one with the smallest H c i(a) (Fig. [I]). At a — 7.5° 
the solutions 9i = 35° and 9 2 = 86° give the same H c \ ((Fig. p, indicating coexistence of 
DVL with vortex lines oriented in these directions. The region 9\ < 9 < # 2 is forbidden for 
vortex-line tilt angles. For moderate anisotropy there is a single 9 and H c i(a) for each a as 
show in Fig. 

ii) Dilute arrangement of vortex-line chains (DVLC). First we calculate 6^(9) by numeri- 
cally minimizing the energy per vortex line of a single chain with respect to the chain period, 
Li, using a simulated annealing algorithm based on a fast convergent series for the energy 

13|jl4}| . The results for e c h(9) and for Li are shown in Fig. |3|. Our results are similar to 



those obtained in Ref. ||. 

It turns out that £ c h(9) is small compared to s s f(9) for all 9: less than 6% for strong 
anisotropy and less than 0.5% for moderate anisotropy. However, as we shall discuss in 
detail shortly, this small correction has non-trivial consequences for strong anisotropy. 

It is interesting to analyze in some detail the chain structure. In Fig. |3] we show the 
interaction energy of a vortex-line pair and the equilibrium positions of the six nearest 
neighbors of a given vortex line for strong anisotropy and 9 = 70°. We note that the chain 
period is considerably smaller than the distance where the vortex-line pair interaction is 
minimum. This results from the long range of the attractive interaction between vortex 
lines. Indeed, in order to reproduce our calculated value for e c h in this case it is necessary 
to add the contributions of all six neighbors shown in Fig. |3]. 

We accurately fit our numerical data for e c h(9) to an analytical function and obtain the 
solutions of Eqs. (j^j and (0) with e(9) = £ s f(9) + £ c h(9) following the same steps as in i). 

For moderate anisotropy we find that the 9(a) and H cl (a) for the DVLC are practically 
identical with the same quantities for the DVL. This results from the very small £ c h(#) in 
this case. 

For strong anisotropy, we find that 9(a) for the DVLC is multivalued in a range of 
a values, with three values of 9 and H c \(a) for each a, similarly to the DVL case. The 
smallest H cl (a), and the corresponding 9(a), are the full curves shown in Fig. f|. Comparing 
the results for the two configurations we conclude that for a < 7.9° H c i(a) for DVLC is 
smaller than that for lines, whereas for a > 7.9° there is practically no difference between the 
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H c i(a) for these configurations (the same is true for a ~ 0°). At a = 7.9° there is coexistence 
between a DVLC arrangement with # 3 = 61° and a DVL one with 64 = 9\ = 86°. The region 
9 3 < 9 < #4 is forbidden for vortex-line tilt angles. 

In conclusion then, we determine if a dilute arrangement of vortex-lines or of vortex-line 
chains is the equilibrium vortex configuration at lower critical field by calculating H c \ for each 
case and by identifying the equilibrium configuration with that for which H c \ is the smallest. 
We find that the vortex-line chain configuration is favored for strong anisotropy and small 
external field tilt angles (< 7.9°). For larger tilt angles the vortex lines are nearly parallel to 
the a-b plane, where there is little difference between the two configurations. For moderate 
anisotropy we find no significant difference between the two H c \. The values chosen in our 
calculation of the anisotropy parameter 7 for what we call moderate and strong anisotropy 
are typical of YBCO and BSCCO, respectively. The parameter k for moderate anisotropy is 
typical of YBCO, but for strong anisotropy is about five times smaller than those typical of 
BSCCO 0. This difference does not alter significantly the above stated conclusions because 
the chain interaction energy does not depend on k. Only the self-energy does. 
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FIGURES 




FIG. 1. Dilute vortex-line arrangement for strong anisotropy. a) 9(a) curve. Dashed line: 
solution of Eq. (|8|). Arrow indicates coexistence point (a = 7.5°). Full line: physical solution 
with smallest H c \. b) Smallest H c \(a) curve. Insets: detail of region where 9(a) and H c \(a) are 
multivalued. Arrows in b) indicate point where H c \(a) curve changes branches. 
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FIG. 2. Dilute vortex-line arrangement for moderate anisotropy. 
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FIG. 3. Numerical results for chains. Triangles: moderate anisotropy. Squares: strong 

anisotropy. a) Chain period. Inset: interaction energy of a vortex-line pair along the x-direction 
(full line) and chain nearest neighbors positions (full circles), b) Chain interaction energy per 
vortex line. 
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FIG. 4. Dilute vortex-line chain and vortex-line arrangements for strong anisotropy compared. 
Dashed lines: (smallest value) dilute vortex-line arrangement (same as in Fig. [|. Full line: 9(a) 
and H c \(a) (smallest value) for dilute vortex-line chain arrangement. Arrow indicates coexistence 
point (a = 7.9°). 
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